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SO.  ABSTRACT  fCanllmM  an  nwh  aM,  II  «•«*••»»  «"*  WMIIIJr  ky  Mack  astronomical  technique 

known  as  speckle  interferometry  shows  considerable  promise  for  satellite  survell-] 
lance  and  imaging.  With  speckle  interferometry,  it  is  possible  to  remove  image 
degradation  caused  by  atmospheric  turbulence.  In  this  way,  resolution  on  satelllt^ 
images  could  be  improved  by  a  factor  of  fifty  with  existing  telescopes.  He  used 
this  method  to  determine  the  mice  and  shape  of  a  high  altitude  satellite.  With¬ 
out  speckle  interferometry,  such  information  would  be  impossible  to  measure. 

•  to  produce  routine  satellite  images  we  are  developing  improved  computer 
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, resolve  high  altitude  satellites.  To  provide  satisfactory  Imagery  of  these 
objects  we  need  telescopes  much  larger  than  existing  instruments.  To -this 
end  we  are  working  to  adapt  a  large  multiple  mirror  telescope  system  for 
speckle  interferometry.  This  telescope  is  the  prototype  for  telescopes,  a 
factor  of  ten  larger  than  any  now  existing. 
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SUmARY 


Turbulence  in  the  Barth's  atiaosphere  causes  "seeing"  which  badly  degrades 
images  of  objects  outside  the  ataosphere.  The  technique  known  as  speckle 
Interferometry  has  shown  considerable  promise  for  recovering  the  information 
lost  to  this  degradation.  We  have  developed  a  television  detector  to  record 
the  short  exposure  data  needed  as  input  for  speckle  interferometry,  and  we 
have  perfected  computer  methods  to  process  this  data.  We  have  performed  a 
test  of  our  system  on  a  satellite  test  target,  successfully  measuring  the 
known  size  and  shape  of  the  object.  Yet  we  still  need  to  Improve  our  ability 
to  track  satellite  targets  so  we  can  observe  larger,  more  interesting  objects 
and  demonstrate  a  complete  image  reconstruction.  Finally,  we  have  pointed  out 
that  telescopes  larger  than  any  existing  are  needed  to  obtain  useful  resolu¬ 
tion  on  high-orbit  objects.  The  new  Multiple  Mirror  Telescope  in  Tucson,  AZ 
will  prove  to  be  important  in  showing  that  a  telescope  as  large  as  25  meters 
can  be  built  and  will  work  to  reconstruct  images  appropriate  to  its  full  2S 
meter  size. 
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I.  INTRODUCTION 


"Seeing, "  caused  by  turbulence  in  the  Earth's  atmosphere,  degrades  images 
from  ground-based  telescopes  to  1  arc-second  resolution  (about  the  size  of  a 
dime  seen  at  one  mile).  Without  this  degradation,  the  resolution  obtainable, 
which  is  said  to  be  "diffraction  limited,”  is  Inversely  proportional  to  the 
telescope  diameter.  For  the  largest  existing  ground  based  telescopes  (*  5 
meter  diameter),  potential  resolution  is  about  0.021  arc-second,  a  factor  of 
50  increase  over  the  seeing  limit.  To  put  this  in  scale,  we  present  in  Table  1 
the  linear  size  of  seeing  and  diffraction  limited  Images  for  a  variety  of 
telescopes  and  satellite  ranges.  Clearly,  much  could  be  gained  by  devising 
means  to  achieve  the  full  resolving  power  of  ground-based  Imaging  systems. 


Table  1.  Telescope  Resolution  Capabilities 


Distance 

(km) 

Atmospheric 

Limit 

1  arc-sec 

4-meter 

Limit 

0.028  arc  sec 

25-meter 

Limit 

0.004  arc  sec 

150 

72  cm 

2  cm 

0.2  cm 

1000 

5  m 

13  cm 

2  cm 

10,000 

50  m 

1.3  m 

20  cm 

40,000  (sync  orbit) 

200  m 

5  m 

SO  cm 

1.5  x  10®  (sun) 

700  km 

20  km 

3  km 

670  kpc  (M31) 

3. 1  pc 

0.089  pc 

0.013  pc 

Astronomers  have  struggled  for  centuries  with  the  desire  to  Improve 
imaging.  Choosing  the  best  sites  and  designing  telescopes  to  minimize  sus¬ 
ceptibility  to  turbulence  can  Improve  images  by  up  to  a  factor  of  two.  In  the 
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1920's,  A.  Michelson  shoved  that  lntarfaroaatry  could  ba  uaad  to  achleva  the 
full  thaoratlcal  raaolutlon  on  astronoalcal  objects.  However,  until  1970,  tha 
hope  of  obtaining  tha  full  raaolving  powar  of  a  large  optical  teleecope  aeeaed 
out  of  reach*  In  1970,  A*  La bey rl a,  a  French  astronomer,  proposed  using  vary 
short  exposure  photographs,  in  a  fora  of  lntarfaroaatry  which  he  called 
"speckle  inter far one try,"  to  realise  tha  full  capabilities  of  existing  ground- 
based  telescopes*  In  1972,  Gesarl,  et  al*.  demonstrated  that  slses  and  shapes 
(but  not  actual  laagas)  could  be  determined  down  to  the  5  aster  (200  inch) 
Faloaar  telescope  diffraction  limit  of  0*021  arc-seconds*  In  1975  Lynds,  et 
al.  showed  that  speckle  lnterferoaetry  could  be  adapted  to  produce  actual 
laagea  down  to  the  0.028  arc-second  Halt  of  the  Kitt  Peak  National 
Observatory  4  aster  (158  inch)  telescope  for  bright  stars. 

In  view  of  its  possible  relevance  to  satellite  identification  and  sur¬ 
veillance,  the  Air  Force  Geophysics  Laboratory  began  a  program  to  apply 
speckle  lnterferoaetry  to  satellites*  It  was  necessary  to  develop  two  areas. 
Because  satellites  are  such  fainter  than  the  brightest  stars,  instrumentation 
to  observe  faint  objects  was  developed.  In  addition,  tha  data  reduction 
methods  were  not  usable  for  general-purpose  image  reconstruction.  We  have  now 
solved  both  of  these  problems*  In  this  report  we  provide  a  non-technlcal 
description  of  the  syatea  and  a  demonstration  of  its  use  for  studying  an 
orbital  satellite,  a  Soviet  Molnlya  coaaunlcatlons  satellite.  We  also  provide 
a  brief  dlacusalon  of  the  need  for  larger  telescopea.  In  paper  II  we  provide 
a  detailed  technical  tutorial  of  current  aethoda  and  status  of  speckle 
lnterferoaetry. 
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II.  HON  SPECKLE  INTERFEROMETRY  WORKS 


Before  we  continue,  It  le  necessary  to  review  the  basic  features  of 
Fourier  Matheaatice  which  are  essential  in  any  discussion  of  optical  system. 

\n  image  is  generally  represented  as  a  set  of  light  Intensities  as  a  function 
of  the  "real”  coordinates  x  and  y.  This  image  may  equally  well  be  represented 
as  a  set  of  harmonic  frequencies,  usually  sines  and  cosines.  In  Figure  1  we 
show  this  graphically*  Assume  a  real  function  f(x)  ■  2x  sin(6*x),  which  is  a 
simple  sine  curve  showing  three  cycles  between  zero  and  one  with  an  amplitude 
of  2.  The  Fourier  representation,  or  transform,  of  this  function  has  zeros  at 
all  points  except  at  a  frequency  of  3  where  (2*f  ■  6ir)  the  value  would  equal 
2,  the  amplitude  of  our  sine  curve.  Yet,  there  is  an  ambiguity  In  this 
representation.  Suppose  we  added  a  constant  in  the  sine  term  with  f(x)  ■  2x 
sln(6*x  +  1.57).  This  function  has  the  same  frequency  representation  as  the 
original  one,  yet  it  is  shifted  relative  to  it.  To  represent  this  in  the 
Fourier  transform,  each  frequency  has  an  amplitude  as  described  above,  and  a 
phase  which  represents  the  shift;  in  this  case  1.57  ■  s/2  •  90s,  so  the  phase 
Is  */2  or  90*.  Since  phase  is  a  circular  function,  it  repeats  at  multiples  of 
2s.  That  is,  a  phase  of  3s  is  the  same  as  a  phase  of  s  etc. 

All  Images  may  be  represented  by  their  Fourier  Transforms.  Most  optical 
systems  act  as  a  form  of  analog  Fourier  Transform  device,  thus  the  usefulness 
of  Fourier  analysis.  Generally  the  amplitude  is  far  easier  to  measure,  in  the 
form  of  a  power  spectrum,  and  the  phase  is  lost.  However,  for  reconstruction 
of  Images  from  Fourier  Transforms,  phase  Information  is  vitally  Important  and 
must  be  recovered  by  some  means.  It  is  intuitive  that  a  large  object  la  re¬ 
presented  by  predominantly  low  frequencies  and  a  small  item  by  high  frequen¬ 
cies.  Thus,  large  objects  tend  to  be  "small"  in  their  Fourier  transforms  and 
vice  versa. 

At  the  image  plane  of  a  telescope  we  measure  the  Fourier  amplitude  of  the 
light  waves  entering  the  telescope  aperture.  In  Figure  2  we  show  this  process 
schematically*  Light  waves  from  a  distant  star  enter  the  telescope  in  Figure  2a 
in  a  form  known  as  a  plane  wave,  which  is  very  much  like  an  ocean  wave  with 
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the  wave  amplitude  and  the  location  of  wave  peaks  and  troughs  the  same  across 
the  whole  aperture*  In  other  words,  the  phase  of  the  Incoming  wave  Is  con¬ 
stant*  The  resulting  Image  is  said  to  be  "diffraction  limited***  Since  the 
image  is  the  Fourier  amplitude  of  the  incoming  wave,  it  Is  easy  to  see  that  a 
larger  aperture  telescope  should  produce  a  smaller  (i*e*,  better  resolution) 
image*  However,  one  gets  perfect  plane  waves  only  in  a  laboratory  or  above  the 
Earth's  atmosphere*  In  practice,  light  waves  passing  through  the  Earth's 
atmosphere  are  broken  both  in  phase  and  amplitude  into  small  segments  of  vary¬ 
ing  extent  averaging  about  10  cm  in  size  as  shown  in  Figure  2b*  Moreover, 
this  breakup  process  changes  very  rapidly;  significant  changes  occur  in  less 
than  1/20  sec  under  typical  conditions*  The  resulting  image,  as  the  Fourier 
transform  of  this  function,  contains  mostly  low  frequencies  which  carry  pri¬ 
marily  information  about  the  breakup  process*  The  image  is  spread  out  con¬ 
siderably  although  it  still  contains  some  high-frequency,  high-resolution 
detail*  Conventional  telescope  photos  have  long  exposures  compared  to  the 
atmospheric  change  time  and  this  remaining  high-frequency,  high-resolution 
detail  Is  washed  out*  Vhat  is  left  is  an  average  image  with  one-arc**second 
resolution  typical  of  a  10  cm  telescope*  Images  are  no  better  for  even  the 
largest  telescopes*  This  image  degrading  process  is  called  "seeing"  by 
astronomers* 

In  1970,  A*  Labeyrie  pointed  out  that  the  small  amplitude,  high-frequency 
information  at  the  telescopes'  full  resolving  power  could  be  recovered  from 
exposure  photos  which  were  shorter  than  the  atmospheric  change  time*  In  these 
short  exposure  photos,  atmospheric  turbulence  is  "frozen"  and  the  telescope 
behaves  like  the  "multiple  aperture  interferometer"  shown  in  Figure  2c*  Some 
of  the  small  patches  distributed  over  the  telescope  aperture  will  be  in  phase 
at  a  given  instant*  Those  in  phase  act  in  concert  to  produce  a  telescope 
diffraction-limited  image,  but  since  many  groups  of  such  concordant  patches 
exist  and  have  different  phases,  a  large  set  of  such  images  is  produced  with 
each  one  shifted  relative  to  another*  Short-exposure  star  photos  in  Figure  3 
show  this  "speckled"  effect*  The  individual  dots  may  approach  perfect  tele¬ 
scope  images:  for  the  binary  star  in  Figure  3c  they  are  double,  and  for  the 
large  star  in  Figure  3a  they  are  larger  than  for  the  unresolved  star  In  3b* 
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This  fact  was  used  In  1975  by  R.  Lynda,  S.  Worden  and  J.  Harvey  at  the  Kitt 
Peak  National  Observatory  to  put  together  the  first  Image  of  the  surface  of 
another  star,  the  supergiant  star  Betelgeuse*  The  resemblence  of  these  photos 
to  laser  speckle  photos  has  led  to  the  process  being  called  "speckle 
interferometry**4  For  complicated  objects  like  near-Earth-orbit  satellites, 
the  speckles  overlap  and  single  Images  are  impossible  to  disentangle*  For 
faint  objects  like  high-orbit  satellites,  only  a  few  photons  are  detectable 
and  individual  speckles  are  not  well-enough  defined  to  be  picked  out* 

Labeyrie  developed  a  method  to  process  this  data  to  produce  the  Fourier 
amplitude  of  the  image*  From  such  data,  size  and  shapes  could  be  determined, 
but  not  actual  images*  For  that,  the  Fourier  phase  is  required  as  well* 
Labeyrlefs  method  was  also  difficult  to  calibrate  to  produce  good  quantitative 
results*  Our  group,  Including  the  authors  of  this  paper,  has  developed  a 
method  to  calibrate  Labeyrie1 s  Fourier  analyses  to  get  accurate  sizes  and 
shapes* 

In  1978,  at  the  Environmental  Research  Institute  of  Michigan,  J*  Fienup 
developed  a  method  to  obtain  the  Fourier  phase  from  an  iterative  mathematical 
computer  method  using  only  the  amplitude  as  an  inupt.  We  have  used  this 
method  with  considerable  promise  to  get  diffraction  limited  images  of  simple 
objects  like  binary  stars,  and  are  now  in  the  process  of  applying  it  to  more 
complex  objects* 
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III.  INSTRUMENT  AND  DATA 


The  photos  shown  In  Figure  3  were  obtained  in  a  very  simple  manner  and 
they  suffer  from  a  number  of  drawbacks.  A  35  on  camera  was  placed  behind  an 
image  intensifler  to  record  the  short  exposure  photos  on  film.  Since  the  image 
intensifier  is  slow  to  return  to  a  ready  condition,  a  shutter  allows  only 
about  i  frame  per  second  so  that  the  Intensifler  will  settle  down  between 
shots.  Since  the  actual  exposure  time  is  about  1/100  second,  99  percent  of 
the  light  is  lost.  Moreover,  the  resulting  photographs  suffer  from  a  variety 
of  photographic  non-uniformities  and  require  a  time-consuming  process  to  digi¬ 
tize  for  computer  processing  (up  to  5  minutes  per  frame).  For  faint  objects 
with  only  a  few  photons  per  frame,  at  least  10,000  frames  are  needed  to  get  a 
usable  result.  With  the  photographic  system  outline  above,  this  means  over  3 
hours  observing  time  and  40  days  data  processing!  To  speed  up  this  process, 
we  replaced  the  film  with  a  television  camera.  The  television  records  a  frame 
every  1/60  of  a  second.  We  run  the  camera  without  a  shutter  and  the  televi¬ 
sion  output  is  directly  computer  processed  in  a  device  known  as  an  "array 
processor."  This  unit  compares  each  frame  with  previous  ones  and  makes  sure 
we  only  record  a  new  photon  when  it  arrives.  The  Fourier  analysis  may  also  be 
done  at  this  time.  With  this  system  we  can  observe,  record,  and  process  data 
to  obtain  our  result  at  the  telescope  all  in  less  than  5  minutes,  to  do  what 
previously  would  have  taken  months. 

As  a  demonstration  of  our  capabilities,  we  used  our  system  on  a  Soviet 
Molniya  communications  satellite  shown  in  Figure  4.  This  object  has  a  12-hour 
orbit,  during  which  time  it  moves  from  an  altitude  of  about  5000  km  to  a  syn¬ 
chronous  distance  of  40,000  km.  With  a  size  of  about  5  m,  at  these  altitudes 
this  object  should  be  marginally  resolvable  with  large  existing  telescopes 
like  the  University  of  Arizona  90  inch  (2.3  meter)  instrument.  Despite  its 
marginal  suitability  for  study,  we  chose  it  since  there  are  few  objects  large 
enough  and  at  a  high  enough  orbit  so  that  relatively  slow  astronomical  tele¬ 
scopes  can  track  them.  For  further  demonstrations  we  plan  to  use  faster 
telescopes  and  larger  objects. 
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For  our  test,  15  minutes  of  data  were  obtained  of  the  object  on  7  July 
1979  beginning  at  0250  Greenwich  Mean  Time.  The  mean  range  of  this  target 
during  the  observation  period  was  21,500  km.  We  also  obtained  a  similar 
amount  of  data  on  a  point  source  (unresolved)  star  for  comparison  purposes. 
About  4  *  10^  photons  were  detected  for  the  target  for  about  15  photons  per 
frame.  We  processed  a  sample  of  the  data  as  described  above  to  produce  the 
Fourier  amplitude  (power  spectrum)  of  both  the  target  and  the  point  source. 

Our  results  are  displayed  in  Figure  5. 

The  results  require  some  words  of  explanation.  We  did  not  try  to  recover 
the  phase  to  produce  an  actual  image  of  this  object  for  two  reasons.  First, 
the  object  is  only  slightly  larger  than  an  unresolved  object  and  it  would  not 
show  much  structure  at  our  resolution.  Second,  the  target  spins  once  per  sec¬ 
ond.  The  several-minute  observation  period  averages  the  spinning,  so  that 
only  the  pattern  of  the  spinning  solar  panels  would  be  visible,  much  as  the 
spinning  blades  of  a  helicopter.  Since  the  target  is  locked  on  the  sun  and  we 
observed  it  at  an  angle,  we  expect  to  see  only  the  oblong  pattern  of  the  spin¬ 
ning  panels.  What  we  have  displayed  is  tlte  "autocorrelation"  of  the  actual 
image.  The  autocorrelation  is  computed  by  Inverse  Fourier  transforming  only 
the  amplitude  or  power  of  the  image's  Fourier  transform  (Ignoring  the  phase), 
but  it  can  be  better  visualized  by  an  analogy.  A  photo  of  an  object  is  taken 
and  a  copy  made  of  it;  the  original  is  then  laid  on  top  of  the  copy  as  shown 
in  Figure  6.  A  light  is  shone  through  the  two  photos  and  the  amount  transmit¬ 
ted  measured  by  a  photocell.  If  we  measure  the  change  of  this  light  level  as 
we  slide  the  copy  in  one  direction  relative  to  the  original,  we  are  measuring 
the  autocorrelation  along  that  direction.  This  function  is  very  useful  in 
studying  the  size  and  shape  of  an  object.  It  follows  that  the  autocorrelation 
of  an  oblong  object  would  be  oblong  and  that  of  a  round  object  would  fall  to 
zero  outside  the  diameter  of  the  object.  In  this  way,  we  can  use  the  autocor¬ 
relation  to  measure  the  size  and  shape  of  the  test  target.  Since  the  target 
is  only  slightly  larger  than  the  telescope  resolution  limit  (see  Figure  7), 
its  autocorrelation  is  only  slightly  larger  than  the  unresolved  star.  None¬ 
theless,  the  expected  oblong  shape  is  readily  apparent.  From  this  data  we 

+2 

were 'able  to  derive  a  size  of  4  j  meters  for  the  target,  In  excellent 
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agreement  with  Its  actual  size  of  A. 9  meters.  The  uncertainty  In  the  results 
Is  set  by  our  ability  to  remove  the  effects  of  the  telescope  Itself,  for  which 
we  use  the  observations  of  the  point-source  star,  and  which  In  turn  is  limited 
by  the  size  of  the  telescope. 

From  these  highly  promising  results,  we  see  two  directions  for  future 
work.  As  mentioned  previously,  faster  tracking  is  needed  to  observe  lower 
orbits  and  larger  objects.  High-rate  tracking  capability  is  available  on 
several  existing  telescopes,  such  as  a  48-inch  telescope  In  Cloudcroft,  NM. 
Using  these  telescopes  we  plan  demonstrations  of  our  method.  Moreover,  we 
recognize  that  actual  reconstructed  images,  rather  than  autocorrelations,  are 
essential.  Based  on  our  recent  astronomical  results,  we  are  confident  that 
the  Flenup  computer  method  will  give  us  the  phase  information  we  need  to 
completely  recover  images. 

If  we  review  Table  1,  it  is  clear  that  the  largest  existing  telescopes 
(~  5  meters  diameter)  do  not  produce  usable  images  for  synchronous  objects, 
even  with  speckle  interferometry.  A  much  larger  telescope  is  clearly  re¬ 
quired,  perhaps  near  25  meters  in  diameter.  The  University  of  Arizona,  in 
concert  with  the  Center  for  Astrophysics  in  Cambridge,  MA,  has  recently  com¬ 
pleted  an  instrument  which  may  prove  that  a  25-meter  telescope  could  be  con¬ 
structed.  This  telescope,  shown  in  Figure  8,  consists  of  six  separate,  72- 
inch  mirrors,  but  it  can  be  used  for  speckle  interferometry  like  a  single  7- 
meter  telescope.  He  are  currently  involved  in  projects  to  evaluate  this 
"Multiple  Mirror  telescope"  for  image  reconstruction  and  to  study  what  would 
be  entailed  in  scaling  it  up  to  25  meters. 
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Figure  1*  Schematic  Representation  of  the  Fourier  Transform 


In  2a,  an  Ideal  situation  with  an  Incoming  plane  wave  from 
a  point  source  Is  shown.  In  2b,  the  atmosphere  has  broken 
the  plane  wave  Into  10  cm  segments  which  may  be  modeled 
Instantaneously  as  a  multiple  aperture  Interferometer  (2c) . 


Figure  2.  The  Image-Forming  Process 
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Figure  3.  Short  Exposure  Speckle  Photos  for  Betelgeuse 
Point  Source,  and  Binary  Star 
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HOTS:  Before  extensive  digital  processing,  autocorrelations  tiers 
computed  in  a  aanner  slailar  to  that  shown  here. 


Figure  6.  An  Iaege  Autocorrelation 
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NOTE:  This  telescope  may  be  used  with  speckle  interferometry  to 
provide  the  resolution  of  a  7  meter  telescope. 

Figure  8.  Multiple  Mirror  Telescope,  Mt.  Hopkins,  Arizona 
(recently  completed) 


